Residential segregation and cultural dissemination: An Axelrod-Schelling model. 
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In the Axelrod's model of cultural dissemination, we consider mobility of cultural agents through the intro- 
duction of a density of empty sites and the possibility that agents in a dissimilar neighborhood can move to them 
if their mean cultural similarity with the neighborhood is below some threshold. While for low values of the 
density of empty sites the mobility enhances the convergence to a global culture, for high enough values of it the 
dynamics can lead to the coexistence of disconnected domains of different cultures. In this regime, the increase 
of initial cultural diversity paradoxically increases the convergence to a dominant culture. Further increase of 
diversity leads to fragmentation of the dominant culture into domains, forever changing in shape and number, 
as an effect of the never ending eroding activity of cultural minorities. 

PACS numbers: 87.23.Ge, 89.20.-a, 89.75.Fb 



The use of agent-based models (ABM) j 1 1 in the study of 
social phenomena provides useful insights about the funda- 
mental causal mechanisms at work in social systems. The 
large-scale (macroscopic) effects of simple forms of (micro- 
scopic) social interaction are very often surprising and gener- 
ally hard to anticipate, as vividly demonstrated by one of the 
earliest examples of ABM, the Schelling J2, Ht] model of ur- 
ban segregation, that shows how residential segregation can 
emerge from individual choices, even if people have fairly 
tolerant preferences regarding the share of like persons in a 
residential neighborhood. 

To gain insights on the question of why cultural differences 
between individuals and groups persist despite tendencies to 
become more alike as a consequence of social interactions, 
Axelrod |j4J proposed an ABM for the dissemination of cul- 
ture, that has subsequently played a prominent role in the in- 
vestigation of cultural dynamics. Questions concerning the 
establishment, spread and sustainability of cultures, as well 
as on the "pros and cons" of cultural globalization versus the 
preservation and coexistence of cultural diversity, are of cen- 
tral importance both from a fundamental and practical point 
of view in today's world. 

The Axelrod model implements the idea that social influ- 
ence is "homophilic", i.e. the likelihood that a cultural fea- 
ture will spread from an individual to another depends on 
how many other features they may have already in common 
10]. The resulting dynamics converges to a global monocul- 
tural macroscopic state when initial cultural diversity is below 
a critical value, while above it homophilic social influence is 
unable to inforce cultural homogeneity, and multicultural pat- 
terns persist asymptotically. This change of behavior has been 
characterized JH H 0, [H as a non-equilibrium phase tran- 
sition. Subsequent studies have analyzed the effects on this 
transition of different lattice or network structures fl fioll . the 
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presence of different types of noise ("cultural drift") HEE!, 
as well as the consideration of external fields (influential me- 
dia) fl3ll and global or local non -uniform couplings lfl4ll . Up 
to now, no investigation of the effects of agent mobility on 
cultural transmission has been carried out, with the exception 
of 11511 . where individuals move following the gradient of a 
"sugar" landscape (that they consume) and interact culturally 
with agents in their neighborhood, i.e., mobility is not cultur- 
ally driven. 

In this paper we incorporate into the Axelrod cultural dy- 
namics the possibility that agents living in a culturally dissim- 
ilar environment can move to other available places, much in 
the spirit of the Schelling model of segregation. This requires 
the introduction of a density of empty sites h in the discrete 
space (lattice) where agents live. As anticipated by lfl5ll the 
expectations are that the agents mobility should enhance the 
convergence to cultural globalization, in the extent that it acts 
as a sort of global coupling between agents. It turns out that 
these expectations are clearly confirmed when the density h of 
empty sites is low enough so that the set of occupied sites per- 
colates the lattice: The transition value depends linearly with 
the number of agents, so that in an infinite system (thermody- 
namical limit) only global cultural states are possible. How- 
ever, for large enough values of h, new phenomena appear 
associated to the mixed Axelrod-Schelling social dynamics, 
including a new multicultural fragmented phase at very low 
values of the initial cultural diversity, a (seemingly first or- 
der) transition to cultural globalization triggered by mobility, 
and the fragmentation of the dominant culture into separated 
domains that change continuously as the result of erosive pro- 
cesses caused by the mobility of cultural minorities. 

In the Axelrod model of cultural dissemination, a culture 
is modelled as a vector of F integer variables {07} (/ = 
1, ...,F), called cultural features, that can assume q values, 
Of = 0, 1, ...q — 1, the possible traits allowed per feature. 
At each elementary dynamical step, the culture {07(1)} of an 
individual i randomly chosen is allowed to change (social in- 
fluence) by imitation of an uncommon feature's trait of a ran- 
domly chosen neighbor j, with a probability proportional to 
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the cultural overlap w.y between both agents, defined as the 
proportion of shared cultural features, 

/=i 

where 6 XtV stands for the Kronecker's delta which is 1 if x = y 
and otherwise. Note that in the Axelrod dynamics the mean 
cultural overlap u)j of an agent i with its ki neighbors, defined 

as 

= ^ >J w y i ( 2 ) 

J=l 

not always increases after an interaction takes place with a 
neighboring agent: indeed, it will decrease if the feature 
whose trait has been changed was previously shared with at 
least two other neighbors. 

To incorporate the mobility of cultural agents into the Axel- 
rod model, two new parameters are introduced, say the density 
of empty sites h, and a threshold T (0 < T < 1), that can be 
called intolerance. After each elementary step of the Axel- 
rod dynamics, we perform the following action: If imitation 
has not occurred and uj^ ^ 1, we compute the mean overlap 
(O and if tDj < T, then the agent % moves to an empty site 
that is randomly chosen. Finally, in the event that the agent i 
randomly chosen is isolated (only empty sites in its neighbor- 
hood), then it moves directly to an empty site. 

We define the mobility m, of an agent i as the probability 
that it moves in one elementary dynamical step (provided it 
has been chosen): 

rrn = (1 - Qi) 0(T - u>i) , (3) 

where Q(x) is the Heaviside step function, that takes the value 
1 if x > 0, and if x < 0. For an isolated agent, that moves 
with certainty, one may convene that its mean cultural overlap 
is zero, so that expression ^ applies as well. The average 
mobility m of a configuration is the average of the mobility of 
the agents: 
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where N is the total number of cultural agents. We will con- 
sider below two-dimensional square lattices of linear size L, 
so that JV = (1 — h)L 2 , periodic boundary conditions, and von 
Neumann neighborhoods, so that the number ki of neighbors 
of an agent i is < ki < 4. We fix the number of cultural 
features to F = 10, and vary the parameters q, h and T, as 
well as the linear size L of the lattice. 

For the initial conditions for the cultural dynamics, N cul- 
tural agents are randomly distributed in the L x L sites of 
the square lattice, and randomly assigned a culture. The sim- 
ulation is stopped when the number n a of active links (i.e., 
links such that < cjy < 1) vanishes. The results shown be- 
low are obtained by averaging over a large number (typically 
5 • 10 2 — 10 4 ) different initial conditions. 




1 2 3 4 5 6 



q/N 

FIG. 1: Order parameter (S m ax}/N versus scaled initial cultural di- 
versity q/N for a very small density of empty sites h — 0.05 and 
different values of the intolerance T — 0.3, 0.7, and of the lattice 
linear size L = 20, 30, 40, as indicated in the inset. 



The usual order parameter for the Axelrod model is 
(Smax) /N, where (S max ) is the average number of agents of 
the dominant (most abundant) culture. Large values (close 
to unity) of the order parameter are the signature of cultural 
globalization. In Fig. Q] we plot the order parameter versus 
the initial cultural diversity scaled to the population size, q/N, 
for a small value of the density of empty sites h = 0.05, and 
different values of the intolerance T and of the linear size L. 
We observe the collapse in a single curve of the graphs cor- 
responding to different lattice sizes and, moreover, that the 
results are rather insensitive to the intolerance values. 

For a fixed value of the initial cultural diversity q, the larger 
the size N of the population is, the more likely an agent can 
share a cultural feature with someone else in the population. 
Hence, as mobility allows contacts with virtually anybody, the 
increase of the population size enhances the tendency towards 
cultural globalization, and the monocultural (ordered) phase 
extends up to higher values of the parameter q. The critical 
value q c of the transition between consensus and a disordered 
multicultural phase diverges with the system size q c ~ N, so 
that in the thermodynamical limit only global cultural states 
are possible for a small density h of empty sites. 

We will focuss hereafter on larger values of the density h 
of empty sites, a regime where the cultural dynamics shows 
strikingly different features. At very low values of the initial 
cultural diversity q (so that cultural convergence is strongly fa- 
vored), the asymptotic states are characterized by low values 
of the order parameter (S m ^ x ) /N. The reason for the absence 
of cultural globalization in this regime is the formation of dis- 
connected monocultural domains, a fact that requires values 
of the density 1 — h of cultural agents (at least) close to (or be- 
low) the site percolation threshold value for the square lattice 
(0.593). This new kind of macroscopic multicultural state is 
thus of a very different nature from the multicultural phase of 
the original Axelrod model (h = 0). The values of the order 
parameter in this fragmented phase, represented in Fig. [2^ as 
a function of q/N with h = 0.5 and T = 0.7 and for sev- 
eral values of L, decrease with increasing lattice size, and the 
expectation is that the order parameter vanishes in the thermo- 
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FIG. 2: (color online) Order parameter {S ma i)/JV versus scaled 
initial cultural diversity q/N for an intermediate value of the den- 
sity of empty sites h = 0.5. Panel (a) corresponds to a high 
value of the intolerance T — 0.7, and different lattice linear sizes 
L = 20, 30, 40, 50, while in panel (b) L = 40, and different val- 
ues of the intolerance T = 0.2, 0.4, 0.7, 0.9 are used. See the text 
for further details. 



dynamical limit, because the largest cluster size below perco- 
lation should be independent of the lattice size. 

The increase of q from the very small values that corre- 
spond to the fragmented multicultural phase has the seemingly 
paradoxical effect of increasing the order parameter (S msa ) /N 
values, i.e., the increase of the initial cultural disorder pro- 
motes cultural globalization. To understand this peculiar be- 
havior, one must consider the effect of the increase of q in 
the initial mobility of the agents. One expects that the higher 
the value of q is, the lower the initial values of the cultural 
overlap u>ij among agents are, and then the higher the initial 
mobility of agents should be. Under conditions of high mo- 
bility, the processes of local cultural convergence are slower 
than the typical time scales for mobility, so that the agents 
can easily move before full local consensus can be achieved, 
propagating their common features, and enhancing the social 
influence among different clusters. In other words, the attain- 
ment of different local consensus in disconnected domains is 
much less likely to occur, and one should expect the coarsen- 
ing of a dominant culture domain that reaches a higher size. A 
straightforward prediction of this argument is that one should 



FIG. 3: (color online) Average mobility m versus time t for h = 0.5, 
L — 30, T = 0.7 and different values of the scaled initial culture 
diversity q/N as indicated. Unlike the other figures, in this case each 
curve represents the results of a single realization. See the text for 
further details. 



observe higher values of (S m . dx ) /N for higher values of the in- 
tolerance T, because agents mobility is an increasing function 
of this parameter (see eq. (O). The numerical results shown 
in Fig.[2j) for different values of T and h = 0.5 nicely confirm 
this prediction, in support of the consistency of the previous 
argument. Interestingly, for high values of the intolerance T, 
an almost full degree of cultural globalization is reached, as 
indicated by the values (S msK )/N ~ 1 of the order parame- 
ter. In those final states almost all agents belong to a single 
connected monocultural cluster. On the contrary, for very low 
values of T when mobility is not enhanced, multiculturalism 
prevails for the whole range of q values. 

To characterize the passage from the multicultural frag- 
mented phase to global consensus with increasing initial cul- 
tural diversity, we have computed the histograms of the values 
of S max /N at values of q where the order parameter increases. 
The histograms display the bimodal characteristics of a first- 
order transition. In a fraction of realizations, the transient 
mobility is able to spread social influence among the clusters 
so that global consensus is finally reached. This fraction in- 
creases with q, to the expense of the fraction of realizations 
where fragmented multiculturality is reached. 

Further increase of the initial cultural diversity q enhances 
the likelihood of agents sharing no cultural feature with any- 
body else in the finite population. The presence of these cul- 
turally "alien" agents decreases the value of the order param- 
eter and the increase of their number with q is concomitant 
with the transition to multiculturality in the original Axelrod 
model (as well as here, for finite populations). We see in 
Fig- Eh that the increase of the intolerance parameter T shifts 
this transition to higher values of q/N, in agreement with the 
enhancement of the convergence to globalization that T pro- 
duces via mobility, as discussed above. Each alien agent has, 
at all times, a mobility m* = 1, and the average mobility 
cannot decrease in time to zero value when they appear. In 
other words, the asymptotic states of the cultural dynamics 
are no longer characterized by m = 0. The time evolution of 
the average mobility m for particular realizations at h = 0.5, 
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FIG. 4: (color online) As a quantitative measure of this erosion phe- 
nomenon we plot here the stationary value of the averaged frac- 
tion Uq I Smax of isolated individuals of the dominant culture ver- 
sus q/N, for h — 0.5, T — 0.7, and L = 30. The inset shows an 
illustrative configuration where erosion can take place. 

T = 0.7, L = 30 and different values of q/N is shown in Fig 
[3] The value of q/N beyond which the stationary average mo- 
bility is larger than zero signals the appearance of these alien 
cultural agents. 

In addition, the restless character of the alien agents has 
an important effect on the geometry of the dominant culture, 
namely its erosion. As an illustrative example, let us consider 
the situation represented in the inset of Fig. |4] in which an 
agent i of the dominant culture is placed at the frontier of a 
cluster, having a single neighbor of his kind, and assume that 
an alien agent j has moved recently to one of the empty neigh- 
boring sites of i. When agent i is chosen for an elementary 
dynamical step, there is a probability 1/2 of choosing agent 
j for an imitation trial. As = 0, and then uii = 1/2, 
the agent i will move from there to a randomly chosen empty 
site whenever the intolerance parameter is T > 1/2. We see 
that, for this particular situation, the erosion of the dominant 
culture cluster will occur with probability one half. 

Note that the erosion of the dominant culture cluster does 
not change the size S max of the dominant culture. It simply 



breaks it up into separate domains, some of them consisting 
of single (isolated) individuals. These isolated members of 
the dominant culture will eventually adhere to domains, to be 
at a later time again exposed to erosion, and so on. Therefore 
the shape and number of domains of the dominant culture (as 
well as that of the other ones), fluctuate forever. The number 
Uq of isolated dominant culture agents reaches a stationary 
value that results from the balance between erosive and adhe- 
sive processes. To quantify the strength of the eroding activity 
of cultural minorities we show in Fig. [4] the stationary value 

of the averaged fraction ( s "° ) of isolated individuals of the 
dominant culture versus the scaled initial cultural diversity, 
for h = 0.5, T = 0.7, and L = 30. Soon after the transi- 
tion from the fragmented multicultural phase to globalization 
occurs, erosion increases dramatically, largely contributing to 
the large values of the stationary mobility m that characterize 
the multicultural states in the model here introduced. 

In summary, the introduction of agents mobility through 
this segregation mechanism into the Axelrod cultural dynam- 
ics leads to an enhancement of the convergence to cultural 
globalization for small densities of empty sites, while for 
larger densities a new type of multicultural fragmented phase 
appears at low values of the initial cultural diversity q, fol- 
lowed by a new transition to globalization for increasing val- 
ues of q that is triggered by the increase in the initial mobil- 
ity. Moreover, in the genuine Axelrod transition from global 
consensus to polarization, the shape and number of cultural 
domains are here dynamically fluctuating by the competitive 
balance of erosive and adhesive processes associated to the 
agents mobility. 
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